Enhancement of the dielectron continuum in Au+Au collisions at y / s/v r /v=200 GeV 
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The PHENIX experiment has measured the dielectron continuum in y / 3jv"iv=200 GeV Au+Au 
collisions. In minimum bias collisions the dielectron yield in the mass range between 150 and 750 
MeV/c 2 is enhanced by a factor of 3.4±0.2(stat.)±1.3(syst.)±0.7(model) compared to the expecta- 
tion from our model of hadron decays. The integrated yield increases faster with the centrality of 
the collisions than the number of participating nucleons, suggesting emission from scattering pro- 
cesses in the hot and dense medium. The continuum yield between the masses of the cj> and the J/tp 
mesons is consistent with expectations from correlated cc production, though other mechanisms are 
not ruled out. 
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Electron-positron pairs, or dileptons in general, have 
proven to be an excellent tool to study collisions of heavy 
ions at ultra-relativistic energies. Because leptons do 
not interact strongly, emission of dileptons from the hot 
matter created at RHIC should leave an imprint on the 
observed dilepton distributions. Emission from the hot 
matter may include thermal radiation and in-medium de- 
cays of mesons with short lifetimes, like the p meson, 
while their spectral functions may be strongly modified. 
However, below the mass of the <f> meson, these sources 
compete with a large contribution of e + e _ -pairs from 
Dalitz decays of pseudoscalar mesons (7r° , 77, 77') and de- 
cays of vector mesons (p, u), 4>). Above the cj> meson mass 
up to 4.5 GeV/c 2 , competing sources are dilepton decays 
of charmonia (J/ip, i/j') and semileptonic decays of D and 
D mesons, correlated through flavor conservation, which 
lead to a continuum of masses. In addition to thermal 
radiation, energy loss of charm quarks in the medium 
might modify the continuum yield in this mass region. 

The discovery of a large enhancement of the dilepton 
yield at masses below the </> meson mass in ion-ion colli- 
sions at the CERN SPS [l[ has triggered a broad theoreti- 
cal investigation of modifications of properties of hadrons 
in a dense medium and of how these modifications re- 
late to chiral symmetry restoration 0] ■ These theoretical 
studies will benefit from the availability of more precise 
data from CERN 0,11] and GSI [5]. An enhanced yield 
was also observed at higher masses, above the (j> meson 
mass Q. Recent NA60 data suggest that the enhance- 
ment can not be attributed to decays of D-mesons but 
may result from prompt production, as expected for ther- 
mal radiation Q. 

The PHENIX experiment at the Relativistic Heavy Ion 
Collider (RHIC) extends these measurements in a new 
energy regime by exploring Au+Au collisions at a cen- 
ter of mass energy of v / saT/v=200 GeV. In this paper we 
present results from minimum bias data taken in 2004. 
Collisions were triggered and selected by centrality using 
beam-beam counters (BBC) and zero degree calorime- 
ters (ZDC). We analyzed a sample of 8xl0 8 minimum 
bias events. 

Electrons and positrons are reconstructed in the two 
central arm spectrometers of PHENIX [8] using Drift 
Chambers (DC), located outside an axial magnetic field, 
which measure their momenta with an accuracy of 
a p /p = 0.7%® l%p/(GeV/c). They are identified by hits 
in the Ring Imaging Cherenkov detector (RICH) and by 
matching the momentum with the energy measured in an 
electromagnetic calorimeter (EMCal) [9J. Electrons are 
reconstructed with an efficiency of ~90%, while a hadron 
contamination of ~20% remains. 

Each central arm covers |At;| < 0.35 in pseudorapidity 
and 7r/2 in azimuthal angle. Because charged particles 
are deflected in the azimuthal direction by the magnetic 
field, the acceptance depends on the momentum and the 
charge of the particle, and also on the radial location of 



the detector component (DC, EMCal and RICH). The 
acceptance for a track with charge q, transverse momen- 
tum pt and azimuthal emission angle <p carL be described 
by: 

1 ^ , . kDC.RICH . , m 

<?>min < <P + q < 9>max (1) 

PT 

where /cdc and fcrticH represent the effective azimuthal 
bend to DC and RICH {k DC = 0.206 rad GeV/c and 
kmcH — 0.309 rad GeV/c). One arm has </> m i n = j^tt 
and max = j|7r, the other arm 4> min = ^ir and 
4>max = j§tt- Only electrons with px > 200 MeV/c are 
used in the analysis. The photon conversion probability 
was minimized by installing a helium bag between the 
beam pipe and the DC, reducing the material to ~0.4% 
of a radiation length. 

In an event the source of any electron or positron is 
unknown and therefore all electrons and positrons are 
combined to pairs, like-sign and unlike-sign. This results 
in a large combinatorial background which must be re- 
moved. The background is computed with a mixed event 
technique, which combines tracks from different events 
that have similar topology (centrality, collision vertex, 
reaction plane). 

In order to achieve the necessary accuracy, all unphysi- 
cal correlations that arise from overlapping tracks or hits 
in the detectors, mostly in the RICH, must be eliminated, 
because they can not be reproduced by mixed events. If 
hits of both tracks of a pair overlap in any detector, the 
event is rejected. About 4% of all pairs are removed by 
this event rejection. Comparing measured like-sign pairs 
with the mixed combinatorial background shows that the 
mixing technique reproduces the shape within the statis- 
tical accuracy of the data. 

The absolute normalization of the unlike-sign combi- 
natorial background is given by the geometrical mean 
of the observed positive and negative like-sign pairs 
2y/N-- N ++ , where, in principle, iV__ and N ++ are the 
measured number of like-sign pairs. There is a small cor- 
related signal also in the observed like-sign pairs, which 
can occur if there arc two c + c~ -pairs in the final state 
of a meson, e.g. double Dalitz decays, Dalitz decays fol- 
lowed by a conversion of the decay photon or two photon 
decays followed by conversion of both photons. These 
"cross" pairs have small masses, typically less than the 77 
mass (550 MeV/c 2 ). 

We therefore determine N and N++ by integrating 

the mixed event distributions after they were normal- 
ized to the 7.5 x 10 6 like-sign pairs measured above 700 
MeV/c 2 . N and N ++ are determined with an accu- 
racy of 0.12%. The normalization is multiplied by 1.004 
to account for the fact that the event rejection removes 
10% more like-sign than unlike-sign pairs. This correc- 
tion factor was estimated, using mixed events, with an 
accuracy better than 50%. Adding the statistical error 
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and the uncertainty due to the event rejection in quadra- 
ture gives an accuracy of 0.25% on the normalization. 

After subtraction of the combinatorial background, 
physical background from photon conversions and cross 
pairs is removed. Because the tracking assumes that the 
e + e~-pair originates at the collision vertex, pairs from 
photons that convert in or outside of the beampipe are 
reconstructed with finite mass and opening angle, which 
is oriented perpendicular to the magnetic field. A cut on 
the orientation of the opening angle in the field removes 
more than 98% of the conversion pairs. 

Cross pairs occur as like and unlike-sign pairs. Monte 
Carlo simulations show that the rate of unlike-sign cross 
pairs accepted in PHENIX is 44% of the rate for like- 
sign cross pairs. To dctcrimc the rate of unlike-sign cross 
pairs, we scale the simulated like-sign cross pair distri- 
bution to the observed like-sign signal, obtained by sub- 
traction of the mixed event background normalized above 
700 MeV/c 2 . We note that the like-sign signal is well de- 
scribed by the Monte Carlo simulation. The simulated 
unlike-sign cross pair distribution is scaled by the same 
factor and subtracted from the unlike-sign signal. The 
uncertainty of this subtraction depends on mass, but is 
< 9% of the final yield. 

Figure [T] shows the mass distribution of c + e~-pairs, 
the normalized mixed event background (B), and the 
signal yield (S) obtained by subtracting the mixed 
event background, the cross pairs and the conversion 
pairs. The insert shows the signal-to-background ratio 
(S/B). The systematic errors (boxes) reflect the er- 
ror on the background subtraction, which is given by 
5s /S = 0.25% • B/S, added in quadrature to the 
uncertainty due to the cross pair subtraction, assumed 
to be 9%S below 600 MeV/c 2 . Despite the small S/B 
ratio, the vector meson resonances u>, <p and J/tp which 
decay directly to e + e~, and an e + c~-pair continuum is 
visible up to 4.5 GeV/c 2 . 

In order to check the background subtraction, a subset 
of data (5xl0 7 events), taken with additional material 
wrapped around the beam pipe to increase the number 
of photon conversions was analyzed. In this data set 
the combinatorial background and the cross pair contri- 
bution is larger by a factor of ~2.5. As shown in Fig. [T] 
the results from both data sets agree well within sta- 
tistical errors, which are 30% in the range from 150 to 
750 MeV/c 2 and much less below. Considering the de- 
creased S/B ratio for the data with the converter we can 
estimate a 0.1% scale uncertainty of the background nor- 
malization, well within the 0.25% systematic uncertainty 
assigned. 

The spectra are corrected to represent the invariant 
yield of e + e~-pairs, with both the e + and e~ in the 
detector acceptance as specified in Eq. [TJ The correc- 
tion is determined using a GEANT simulation [l(| of the 
PHENIX detector that includes the details of the detec- 
tor response. Simulated e + e _ -pairs are reconstructed 



with the same analysis chain and all cuts applied. The 
correction is determined double differentially in px and 
mass of the e + e~-pair. The reduction of the electron 
reconstruction efficiency (0.92+/-0.03) due to detector 
occupancy is corrected for. Systematic uncertainties on 
the correction can be summarized as: (i) 13.4% on di- 
electron reconstruction, which is twice the uncertainty 
on the electron reconstruction efficiency i, (ii) 6% con- 
version rejection cut, (iii) 5% event rejection and (iv) 3% 
occupancy. These uncertainties are included in the final 
systematic error on the invariant e + e _ -pair yield. 

Figure [2] compares the invariant yeild to the ex- 
pected yield from meson decays and correlated decays 
of charmed mesons. The cocktail of hadron decay 
contributions was estimated using PHENIX data for 
meson production when available. As input distribu- 
tions we use the measured n, r), <f>, J /ij) yield and spectra 



U, 12, 13, 14, fla ]. For other mesons we use the tot 
scaling procedure outlined in The systematic uncer- 
tainties depend on mass and range from 10 to 25%. 

For the continuum below the J/tp the dynamic corre- 
lation of c and c is essential, but unknown. We make 
two assumptions: (i) the correlation is unchanged by the 
medium and equal to what is known from p+p collisions. 
In this case we use PYTHIA [l6[ scaled from the P+p 



equivalent cc cross section of 567+57+193 ^ibarn [17 1 
to minimum bias Au+Au collisions proportional to the 
number of binary collisions (258 + 25) [12]. We note that 
the pt distribution for electrons generated by PYTHIA 
is softer than the spectra measured in p+p data but co- 
incides with those observed in Au+Au As a second 
assumption (ii) there is no dynamical correlation, i.e. the 
direction of c and c quarks are uncorrelated. We keep the 
overall cross section and the px distributions fixed to ex- 
perimental data [§{]. Other contributions from bottom 
and Drell-Yan are expected to be small in the mass re- 
gion below the J/tjj peak. Each e + and e~ must fall in 
the PHENIX acceptance, given by Eq. [TJ 

The data below 150 MeV/c 2 are well described by the 
cocktail of hadronic sources. The vector mesons u>, 4> and 
J/ip are reproduced within the uncertainties. However, 
the yield is substantially enhanced above the expected 
yield in the continuum region from 150 to 750 MeV/c 2 . 
The enhancement in this mass range is a factor of 3.4 
± 0.2(stat.) ± 1.3(syst.) ± 0.7(model), where the first 
error is the statistical error, the second the systematic 
uncertainty of the data, and the last error is an estimate 
of the uncertainty of the expected yield. Above the </> 
meson mass the data seem to be well described by the 
continuum calculation based on PYTHIA. This is some- 
what surprising, since single electron distributions from 
charm show substantial medium modifications and 
thus it is hard to understand how the dynamic correla- 
tion at production of the cc remains unaffected by the 
medium. A complete randomization of that correlation 
(see Fig ^ leads to a much softer mass spectrum and 
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FIG. 1: Uncorrected mass spectra of all e + e~ pairs, mixed events background (B) and signal (S) with statistical (bars) and 
systematic (boxes) uncertainties shown separately. The signal from the runs with additional converter is shown with statistical 
errors only. The insert shows the S/B ratio. The mass range covered by each data point is given by horizontal bars. 




FIG. 2: Invariant e + e~-pair yield compared to the yield from the model of hadron decays. The charmed meson decay 
contribution based on PYTHIA is included in the sum of sources (solid black line). The charm contribution expected if the 
dynamic correlation of c and c is removed is shown separately. Statistical (bars) and systematic (boxes) uncertainties are shown 
separately; the mass range covered by each data point is given by horizontal bars. The systematic uncertainty on the cocktail 
is not shown. 



would leave significant room for other contributions, e.g. 
thermal radiation. 

To shed more light on the continuum yield we have 
studied the centrality dependence of the yield in three 
mass windows, below 100 MeV/c 2 , from 150 to 750 



MeV/c 2 and 1.2 to 2.8 GeV/c 2 . The top panel of Fig. 02 
shows the centrality dependence of the yield in the mass 
region 150-750 MeV/c 2 divided by the number of par- 
ticipating nucleon pairs (./V part /2). For comparison the 
yield below 100 MeV/c 2 , which is dominated by low px 
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z^a) Au+Au at = 200 GeV 
Yield /(N parl /2) 

• DATA (1 50<rrj a <750 MeV/c 2 ) 
I ] COCKTAIL {150<rr^ o <750 MeV/c 2 ; 




DATA (0<iT), e <100 MeV/c ) 
COCKTAIL (0<rr^ B <100 MeV/c 2 ) 
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FIG. 3: Dielectron yield per participating nucleon pairs 
(-/Vpart /2) as function of N par t for two different mass ranges 
compared to the expected yield from the hadron decay model. 
The two lines give ±1<t systematic uncertainty. For the data 
statistical and systematic uncertainties are shown separately. 



pion decays, is shown in the lower panel. For both inter- 
vals the yield is compared to the same yield calculated 
from the hadron cocktail. In the lower mass range the 
yield agrees with the expectation, i.e. is proportional to 
the pion yield. In contrast, in the range from 150 to 750 
MeV/c 2 , the observed yield rises significantly compared 
to the expectation, reaching a factor of 7.7 ± 0.6(stat.) 
± 2.5(syst.) ± 1.5(model) for most central collisions. 
The increase is qualitatively consistent with the conjec- 
ture that an in-medium enhancement of the dielectron 
continuum yield arises from scattering processes like 7T7t 
or qq annihilation, which would result in a yield rising 
faster than proportional to N part . 

We normalize the yield in the mass region 1.2 to 2.8 
GeV/c 2 to the number of binary collisions (Fig. [4j) , which 
is the correct scaling for pairs from charmed meson de- 
cays The normalized yield shows no significant cen- 
trality dependence and is consistent with the expectation 
based on PYTHIA. It is also likely that a scenario where 
the correlation between the c and c is randomized will re- 
quire an additional source, e.g. a contribution from ther- 
mal radiation. This contribution could increase faster 
than linear with -/Vp art and therefore the apparent scaling 
with N co u may be a mere coincidence. We note that this 
coincidence may have been observed in this mass region 
at the CERN SPS where a major prompt component 
has now been suggested by NA60 data 0]. 

In conclusion, measurements of Au+Au collisions at 
77^7=200 GeV in the mass range 150-750 MeV/c 2 
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FIG. 4: Dielectron yield per number of collisions N co u in the 
mass range 1.2 to 2.8 GeV/c 2 as function of -/V par t. Statistical 
and systematic errors are shown separately. Also shown are 
two bands corresponding to the two different estimates of the 
contribution from charmed meson decays. The width of the 
band reflect the uncertainty of the charm cross-section only. 



show a significant enhancement of the dielectron con- 
tinuum and exhibit a clear increase with centrality of 
the collision. The observed yield between (j> and J/ip is 
consistent with the expectation from correlated cc pro- 
duction, but does not exclude other mechanisms. 
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